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Aims. We report near-infrared observations of the supergiant donor to the eclipsing high mass X-ray binary pulsar IGR J 18027-20 16. 
We aim to determine its spectral type and measure its radial velocity curve and hence determine the stellar masses of the components. 
Methods. ESO/VLT observations of the donor utilising the NIR spectrograph ISAAC were obtained in the H and K bands. The multi- 
epoch H band spectra were cross-correlated with RV templates in order to determine a radial solution for the system. 
Results. The spectral type of the donor was confirmed as B0-1 I. The radial velocity curve constructed has a semi-amplitude of 
23.8 ±3.1 km s -1 . Combined with other measured system parameters, a dynamically determined neutron star mass of 1.4 ± 0.2 - 1.6 + 
0.3 M G is found. The mass range of the B0-B1 I donor was 18.6 ± 0.8 - 21.8 ± 2.4 M G . These lower and upper limits were obtained 
under the assumption that the system is viewed edge-on (i = 90° with ft = 0.89) for the lower limit and the donor fills its Roche lobe 
(/J = 1 with i = 73.1°) for the upper limit respectively. 
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1. Introduction 

The precise nature of the fundamental properties of neutron star 
(NS) matter has been an on-going field of research for many 
decades. The physical properties of matter that exist in the ex- 
treme densities and pressures found within a NS can be deter- 
mined by the NS equation of state (Eo S). At this junctur e how- 
ever there exist over 100 proposed EoS (Ka per et al.l20 06). Only 
one can be physically valid. Observational data can assist in re- 
ducing this number by eliminating contending EoS that place 
unrealistic constraints on the mass range of observed NSs. 
Many EoS predicting the presence of exotic hadronic matter 
have been ruled out by the recent mass determination of 1.97 
± 0.04 Mp for the binary m illisecond pulsar (MSP) PSR J1614- 
2230 dDemorest et al.l2010t) . In the case of PSR J1614-2230 it is 
believed that the over-massive NS in this system was created by 
the accretion of matter whilst being spun-up to become a MSP. 
Another possible means of high mass NS production may oc- 
cur where massive progenitors have yielded a high pre-SN core 
mass. High Mass X-ray Binary systems (HMXBs) are potential 
testing grounds for this idea as in the majority of cases the pro- 
genitor of the NS will be more massive than the donor star we 
can observe. 

The only means of performing unambiguous NS mass determi- 
nations in X-ray pulsar binaries is by observing eclipsing bi- 
nary systems in which we can constrain the inclination angle. 
There are presently 10 such systems kno wn with 8 i n whic h 
the NS mass has been calculated (e.g. Mason et al. I (I2TJTT1) : 



* Based on observations carried out at the European Southern 
Observatory under programme ID 085.D-0539(A) 



IMason et ail (l2010l) : lOuaintrell et al.l (120031) 1. In this paper we 
present an orbital solution from observations of the donor within 
the HMXB accretion driven pulsar IGR J18027-2016. The mass 
donor within this system is heavily obscured and reddened. 
Using near infrared spectroscopy of the donor star we have cal- 
culated an orbital solution to measure the mass of each of the 
components within the system. 

IGR J 18027-20 16 was first detected by INTEGRAL in 2003, 
dRevnivtsev et al.l20 04) and was spatially associated with the X- 
ray pulsar SAX J18027. 7-2017 which was discovered serendip- 
itously during observations of the LMXB GX9+1 and found 
to have a pulse period of 139.6 s (lAugello etalJl2003h . IGR 
J18027-2016 has a measure d hydrogen co lumn density, Nh = 
9.1 ± 0.5 x 10 22 cm -2 dWalter et alj|2006l) greatly in excess of 
the line of sight column density of Nh = 1.0 x 10 22 crrT 2 . This 
is indicative of intrinsic absorption occurring within the system, 
with matter from the stellar wind of the dono r forming a dens e 
spherical shell around the NS as it is accreted (iHill et al.l l2005). 
XMM-Newton observations in 2004 narrowed the source posi- 
tion down (with an uncertainty of 4") to the location RA(2000.0 ) 
= 18 A 02 m 42.0 5 and Dec = -20° 17' 18" (I Walter et al) [2006). 
Using this improved position Ma setti et al.l d2008l) were able to 
identify the donor as 2MASS J 18024 194-20 17 172 (J = 12.7, H 
= 11.9,K= 11.5). 

From timing analysis IGR J18027-2016 was fo und to have an or- 
bital period of 4.5696 + 0.0009 days (IHill et alj|2005l) . The orbit 
based upon the sinusoidal modulation of pulse arrival times was 
found to be approximately circular. Using a donor mass-radius 
relation togethe r with approximations of the Roche lobe radius 
iHiil et al.l (120051 ) proposed that the donor is a 09 - B 1 supergiant. 
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Table 1. The phase, radial velocity and telluric standard for each IGR J 18027-20 16 spectrum. 



Mid-point of Observations (UT) 


MJD 


Phase 


Radial velocity / km s 1 


Telluric Std 


2010 July 13.161 


55390.16060 


0.073 


51.8 ± 


9.9 


Hip 088201 


2010 July 09.017 


55386.01693 


0.166 


76.9 ± 


9.9 


Hip 089744 


2010 April 22.336 


55308.34560 


0.169 


68.4 ± 


9.9 


RV std 


2010 June 07.210 


55354.21010 


0.205 


85.2 ± 


9.9 


Hip 092931 


2010 Aug 10.197 


55418.19702 


0.208 


79.9 ± 


9.9 


Hip 092322 


2010 Sept 06.995 


55445.99528 


0.291 


76.0 ± 


9.9 


Hip 094859 


2010 June 12.252 


55359.25220 


0.309 


69.2 ± 


9.9 


RV std 


2010 June 22.074 


55369.07022 


0.457 


64.8 ± 


9.9 


RVstd 


2010 Aug 25.171 


55433.17163 


0.485 


47.7 ± 


9.9 


Hip 100170 


2010 Aug 16.111 


55424.11101 


0.502 


48.6 ± 


9.9 


Hip 092470 


2010 Sept 08.080 


55447.07987 


0.529 


33.4 ± 


9.9 


Hip 090804 


2010 July 29.211 


55406.21077 


0.585 


42.6 ± 


9.9 


Hip 089384 


2010 June 09.095 


55356.09537 


0.618 


37.8 ± 


9.9 


Hip 083535 


2010 May 26.392 


55342.39190 


0.619 


22.7 ± 


9.9 


RV std 


2010 Aug 08.104 


55416.10351 


0.750 


19.9 + 


9.9 


Hip 101552 


2010 Aug 31.054 


55439.05354 


0.772 


24.2 + 


9.9 


Hip 094378 


2010 June 19.167 


55366.16684 


0.822 


56.4 ± 


9.9 


HD 169101 


2010 July 08.191 


55385.19053 


0.985 


56.3 ± 


9.9 


Hip 090978 



O ptical and NIR spe ctral analysis of the mass donor performed 
bv lChatv et al.l (120081) detected prominent Paschen and Brackett 
series hydrogen lines in addition to He i and He n lines in 
emission. Combining their spectroscopy with SED modelling 
they proposed that the mass donor was an early B super- 
giant. Obtaining hig her resolution spectra in the I and K bands 
iTorreion et alJ (1201 Oh narrowed the spectral type down to B 1 lab 
- Bl lb. The short orbital period of IGR J 18027-206 of only 4.6 
d eliminated the donor being a B 1 lab supergiant, they proposed 
a classification of B 1 lb for the mass donor lying at a distance of 
12.4 kpc. 



2. Observations and Data Reduction 

Observations were conducted between 2010 May 26th and 2010 
Sept. 08th using the NIR spectrograph ISAAC on the VLT in 
the SW MRes mode utilising a 0.8" slit width. Although the 
mass donor in IGR J18027-2016 is relatively faint (H = 11.91) 
we were able to obtain medium resolution (R ~ 3000) and high 
S/N spectra in the H and K s bands. For the H band multi-epoch 
science and RV template exposures centred on 1 .7 fim were ob- 
tained of 2400s and 200s respectively. For the K s band two sci- 
ence exposures centred on 2.06 and 2.15 //m were taken for a to- 
tal of 2240s. Reduction was performed using the ISAAC pipeline 
with OH skylines used to wavelength calibrate the spectra. The 
resulting data has a count rate less than 10 000 ADU and thus 
no correction for detector non-linearity was necessary. The RV 
template observed was Hip 89262, with a known RV of -8 km 
s _1 , close in spectral type to the target (B0. 5 la) and bright (H ~ 
6.8). 

We obtained a dataset containing 18 spectra. Telluric correction 
was employed to remove atmospheric features from each spec- 
trum. Unfortunately, no telluric standards were obtained for the 
first four spectra in our dataset (see Table 1). In this case we em- 
ployed the RV template (which had an airmass close to that of 
the target) as a telluric standard. In order not to contaminate our 
target spectrum we first removed any non-telluric spectral lines 
from the RV template before telluric correction was attempted. 
The continuum normalised H band spectra comprising our full 
dataset ordered by phase are shown in Fig. 1. 



3. Data Analysis 

Using the standard IRAH3 routine fxcor, radial velocities were 
determined by cross-correlating the region around the He i 1 .70 
/mi absorption line in the 8 science spectra against the same re- 
gion in the RV template spectra which were obtained shortly af- 
ter each science exposure. The resulting radial velocities were 
then corrected to the solar system barycentre and are reported 
in Table 1 . The ecce ntricity of the system is reported as e < 0.2 
(lAugello et al. 1 120031) and al though no precise measurement of 
the eccentricity is reported, Hilletal. (2005) found sinusoidal 
modulation from an analysis of pulse arrival times that indicates 
the NS is in a circular orbit around the mass donor. We have thus 
fitted the radial velocities of the super giant dono r with a sinu- 
soidal solution, using the ephemeris of iHill et all (120051) which 
specifies the epoch of mid-eclipse as 



T(MJD) = 52168.26(4) + 4.5696(9)JV 



(1) 



the uncertainties in brackets refer to the last decimal place 
quoted and N is the cycle number. At the epoch of our observa- 
tions the accumulated uncertainty in phase is 0. 1 1 . We therefore 
fitted our data with two models, one in which the zero phase is 
fixed and the other in which we allowed the zero phase to vary 
as a free parameter. 

In our first model we fixed the zero phase and allowed the 
systemic velocity and RV semi-amplitude to vary as free param- 
eters. Fitting a sinusoid to our data we found the radial velocity 
semi-amplitude Ko = 24.4 + 3.2 km s _1 and systemic velocity 
y = 51.7 ± 2.4 km s _1 . To obtain this solution the uncertainties 
in RV for each data point had to be scaled to + 9.9 km s _1 thus 
reducing chi-squared to unity. 

Alternatively fitting a sinusoid to our data with three free pa- 
rameters (allowing in this case the zero phase to vary) we found 
the radial velocity semi-amplitude Ko = 23.8 ± 3.1 km s _1 , sys- 
temic velocity y = 52.5 + 2.4 km s and a phase shift of -0.03 
+ 0.02. Here the uncertainties were scaled to + 9.7 km s _1 . The 
best-fit phase offset found is within the accumulated phase un- 
certainty of the ephemeris. We prefer this fit and use the data 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Continuum normalised H-band spectra centred on He i 
1.70 fim of IGR J18027-2016 in order of phase from bottom to 
top. 



from it for our subsequent calculations of the system parameters 
of IGR J18027-2016. Both fits are shown in Fig 2. 

The projected semi-major axis of the neutron s tar's orbit 
from X-ray pulse timing delays was measured by iHill et all 
(2005) as a x sin / = 68.0+1 light seconds. From this, the semi- 
amplitude of the neutron star's radial velocity may be calculated 
using 

a x sini = K x (2) 

for a circular orbit to give Kx = 324.4 + 4.7 km s" 1 . 

To determine the component masses of the system we must 
first consider the mass ratio of the system q which is equal to the 
ratio of the semi-amplitudes of the radial velocities for each star 
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Mx^Ko 
Mq K x 



(3) 



where Mx and Mo are the masses of the neutron star and su- 
pergiant star respectively, and Kx and Kq are the corresponding 




Fig. 2. Radial velocity data for the supergiant mass donor in IGR 
Jl 8027-2016. The solid line is the best fitting sinusoid with three 
free parameters, the dashed line is that with a fixed zero phase 
in line with the publis hed ephemeris. T he orbital phase is based 
upon the ephemeris of Hi ll et al.l d2005l) . 



Table 2. System parameters for IGR J18027-2016. 



Parameter 


Value 


Ref. 


Observed 






a x sin i / It sec 


68 ± 1 


[1] 


P/d 


4.5696 ± 0.0009 


[1] 


r 90 / MJD 


52168.26 + 0.04 


[1] 


e 


^0.2 


[2] 


e e /deg 


34.9 ±4.6 


[1] 


K 1 km s" 1 


23.8 + 3.1 


[3] 


Assumed 






SI 


= 1 




Inferred 






K x 1 km s _1 


324.4 ± 4.7 




q 


0.073 ± 0.01 





:6.3 



1.000 
73.1 i 
1.58 ± 0.27 
21.8 ± 2.4 
33.1 ± 1.2 
19.8 + 0.8 

19.8 + 0.7 

[1] Hill et al. (2005); [2] Aueello et al. (2003) 
[3] this paper 



semi-amplitudes of their radial velocities. In addition, for circu- 
lar orbits, 



i'/deg 

Mx I M 
M /M G 
«/Ro 
fi L /Ro 

^o/Ro 



0.893 ± 0.078 
90.0 

1.36 + 0.21 
18.6 + 0.8 
31.4 + 0.5 
18.8 + 0.3 
16.8 ± 1.5 



M = 



K X 3 P 



and similarly 



M x 



2nG sin~ i 
2ttG sin j i \ q 



(4) 



(5) 



where i is the inclination to the plane of the sky and P is the 
orbital period. The system inclination can be found from the ge- 
ometric relation 



sin; 



1/2 



COS ft, 



(6) 
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Fig. 3. Left: IGR J18027-2016 continuum normalised H band spectrum. Right: Continuum normalised K s band spectrum of IGR 
J18027-2016. 



where 6 e is the eclipse half-angle, is the Roche lobe radius of 
the supergiant, B is the ratio the supergiant's radius to that of its 
Roche lobe and a is the separation between the centres of mass 
of the two stars. The Roche lobe radius may be approximated by 



— x A + Blogq + C log 2 q 
a 



(7) 



where the constants have been determined by Joss & Rappaport 
(1984) as 

A * 0.398 - 0.026O 2 + 0.004O 3 (8) 



B * -0.264 + 0.052Q 2 - 0.015Q 3 
C * -0.023 - 0.005Q 2 



(9) 
(10) 



D. is the ratio of the spin period of the supergiant to its orbital 
period. We have assumed that the supergiant is close to Roche 
lobe-filling and will likely be rotating synchronously with the 
orbit, so Q = 1, and the ecl ipse half angle is taken to be 6 e = 
34.9° ± 4.6° dHill et al.ll2005l) . 

Equations (2) - (10) allow the masses of the two stars to 
be determined in two limits. First, assuming that the system is 
viewed edge-on (in which case i = 90°) we can find a lower 
limit to the Roche lobe filling factor B and lower limits to the 
stellar masses. Secondly, assuming that the supergiant fills its 
Roche lobe (in which case B — 1) we can find a lower limit to 
the system inclination i and upper limits to the stellar masses. 
The results in each limit are shown in Table 2. The true values 
of the masses of the component stars within the system may lie 
anywhere between these limits for corresponding combinations 
of B and i. In order to determine the uncertainties on each derived 
parameter, we performed a Monte Carlo analysis which propa- 
gated the uncertainties through 10 4 trials, assuming a Guassian 
distribution in each individual uncertainty. 
Between phases <p ~ 0.4 - 0.6 the He i 1 .70 /mi absorption line 
displays asymmetric features for a number of the spectra (Fig 1). 
To investigate the effect these asymmetries have on the final or- 
bital solution, we omitted spectra between this phase range and 
calculated a new set of NS and donor star masses. We found the 
impact on the derived orbital solution of these asymmetries to 
be negligible and have included the full set of measured radial 
velocities for completeness. We are unsure of the exact physical 



cause of these asymmetric features, but there appears to be a cor- 
relation between the airmass of the observation and the degree 
of asymmetry observed. We favour atmospheric effects arising at 
the time of the observation to be the cause; whilst acknowledg- 
ing that they have a negligible effect upon the calculated system 
masses. 

4. IGR J1 8027-201 6 Spectral classification 

The H band spectrum is displayed in Fig. 3(a). Here we see the 
He i 1.700 /mi line is present but there is no detection of He n 
1.692 /mi, which is seen in absorption in early and late O super- 
giants and typically absent in early B supergiants. Further spec- 
tral indicators include the two hydrogen lines, Brackett 10 (1.736 
/mi) and Brackett 11 (1.681 /mi) which dis play relatively wea k 
absorption features indicative of B0 - Bl I (lHanson et al.l l2005). 
Examining the K, band, the absence of lines due to the C iv 
triplet (2.069, 2.078 and 2.083 /mi) combined with the He 1 2.058 
/mi line in emission implies I GR J18027-2016 is not an O type 
supergiant (Ma son et al.ll2.Q09b . Further constraining the spectral 
type, we can see that He i 2.058 /mi is in emission and when 
combined with the observed He i 2.112 /mi line in absorption, 
this is indicative of a B0-B2 I classification. The N m 2.1 15 /mi 
line seen in emi ssion allows us to again converge on a spectral 
type of B0-B1 1 (lHanson et alj|2005b . 

We can also see that Br y displays relatively weak emission, un- 
fortunately we have no spectral coverage of the He n 2.189 /mi 
line region. As the spectral features in the K, band are not highly 
dependent on luminosity, we have difficulty distinguishing be- 
tween luminosity sub-classes. We can not localise the classifica- 
tion to any better than B0 - Bl I. This is broadly in agreement 
with the spectral type found by Torrei on et al.l (HoiO) of Bl lb. 
IGR J 18027-20 16 does not appear to exhibit any secular changes 
in spectral mo r phol ogy, both the K band spectrum obtained by 
iTorreion et al.l d2010l) and the spectrum shown in Fig. 3(b) are 
separated by 5 years but are very similar. 

5. Conclusions 

Within this paper we have presented the analysis and results of 
multi-epoch observations of the eclipsing high mass X-ray bi- 
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nary IGR J 18027-20 16 performed at the ESO/VLT with the NIR 
spectrograph ISAAC. These have enabled us to make the first 
measurements of the dynamical masses of both system compo- 
nents. Constructing an orbital solution from near-IR radial ve- 
locity measurements and the orbital parameters of the system, 
provides a dynamically determined neutron star mass of 1.4 + 
0.2- 1.6 ±0.3 M G . 

The mass and radius of the supergiant donor, M ~ 18-22 M 
and R ~ 17 - 20 R is in agreement with t hat suggested fro m 
X-ray observations (M ~ 21 M Q , R ~ 19 R^ lHill et alJ (l2005h \ 
We find the spectral type of the donor to be B 0-B1 I, this is 
in acc ord with that found previously (Bl lb) by Torreio n et alj 
(120101) . Although we note that due to the lack of spectral cov- 
erage encompassing the He n 2.189 //m line, this prevented us 
from conducting non-LTE stellar atmosphere modelling and ex- 
tracting the full range of stellar parameters. Particularly use- 
ful would have been a determination of the stellar temperature 
that would enable us to accurately constrain the luminosity sub- 
class. From an examinat ion of a sample of Galactic B super- 
giants dSearle et alj2 008) the donor has a wide possible range of 
luminosities, log(L/L ) ~ 5.5 - 5.7 and temperatures log T e ff ~ 
4.3 - 4. 5. Comparing this parameter range to evolutionary tracks 
(Fig. 8 ISearle et alj (120081) ) leads us to suspect that the donor's 
progenitor mass was ~ 30 M . 
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